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Abstract 

In this letter, we show how one can solve easily the Potts-3 + branching inter- 
actions and Potts-oo matrix models, by the means of the equations of motion (loop 
equations). We give an algebraic equation for the resolvents of these models, and 
their scaling behaviour. This shows that the equations of motion can be a useful tool 
for solving such models. 
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1 Introduction 



Random matrices are useful for a wide range of physical problems. In particular, they 
can be related to two-dimensional quantum gravity coupled to matter fields with a non- 
zero central charge C [Q. While C < 1 models are relatively well understood, the C > 1 
domain remains almost totally unknown : there is a C = 1 "barrier" . When studying C 7^ 
models, we are led to consider multi-matrix models 0, |^ which are often non-trivial. One 
class of difficult matrix models corresponds to the g-state Potts model (in short: Potts-g) 
on a random surface. This model is a g-matrix model where all the matrices are coupled 
to each other, thus making difficult the use of usual techniques such as the saddle point 
or the orthogonal polynomials method. Moreover, the g — > 4 limit corresponds to C ^ 1, 
thus, by solving Potts-g models, we shall gain a new understanding of the C = 1 barrier. 

In this letter, we show that, contrary to what was previously thought, one can use the 
loop equations to solve the Potts-3 random matrix model, and we find that the resolvent 
(which generates many of the operators of the problem) obeys an algebraic equation that 
we write explicitely. 

We also show that this method applies when one adds branching interactions (gluing 
of surfaces, also called "branched polymers") 0] and we derive the critical line of this 
extended model. The extension to the model with branching interactions and the study of 
its phase diagram is necessary to verify [^'s conjecture about the C = 1 transition. 

Finally, we apply the method to the Potts-00 matrix model, which corresponds to 
C = 00. 

As this work was approaching its completion, a paper appeared on the dilute Potts 
model [^], which partially overlaps our present work. In this article, the author also has an 
algebraic equation for the conventional Potts-3 model. Here, we go further as we consider 
the Potts-3 + branching interactions model. Moreover, his method is quite different : while 
he uses analytical considerations on the resolvents and large-A^ techniques, we solve our 
model by the loop equations method, which can be extended to finite problems and is 
also more adapted to the use of renormalization group techniques 0, |^. 



2 The Potts-3 + branching interactions model 

Let us define : 

Z = j ci$e-^'^(*) 
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6 = 6i + S^i. We shall also use later the notation = Id. 



$, 60, 61 and 6-1 are (3A^) x (3A^) and $1, $2, $3 iV x iV hermitean matrices. V is a 
general two-variable function, and will mainly appear through its partial derivatives U and 
c with respect to ^^|^ and ^^l^*^*^ respectively. If these are constants, then we recover 
the conventional Potts-3 model (no branching interactions). This model was given partial 
solution by J.M. Daul |P, by considering the analytic structure of the resolvents. He had 
its critical point and its associated critical exponent. He did not know, however, if the 
resolvent obeyed an algebraic equation. We shall give here the expression of this equation 
for the conventional and extended Potts-3 model. We also derive the critical line of the 
extended model and check it corresponds to Daul's result in the particular case of the 
conventional model. 



Let us note, for convenience : 

th^,■■■i„^^ = ^(tr5i,$5,,$ . . . <5i„$*^) (4) 

where ii, . . ., in can be +1, —1 or 0. This trace is non-zero if and only if ii + ... + in = 
(mod 3). (. . .) is the expectation value of (...) : 

(...) = i/rf$(...)e-^^^(*) (5) 

A trace will be said to be "of degree m" if there are m matrices $ in it. For example, the 
above trace is of degree k + n ~ 1. 

Let us now use the method of the equations of motion (or loop equations). If we make 
the infinitesimal change of variables in Z : 

$_$ + e5,,<l>5,,...$5,„ (6) 

with 

^1 + ^2 + • • • + ^n. = (mod 3) 
then we obtain the expression of the general equations of motion : 

n-l 

g tii...i„<S>'^ + U ti^,„i^ip + C (t(ii+i)i2...(i„-l)$ + t(ii-l)i2...{i„+l)<i>) — X! ^h-ij ^ij+i-in = i'^) 

i=i 

The first three terms come from the transformation of V{^), and the last one, from the 
jacobian of the transformation. 

Eq. d^) relates any expectation value of trace containing a quadratic term (i.e. a 
term) to expectation values of traces of lower degrees. The problem is that we do not 
have any recursion relation for more general expectation values like tii...j„$ where all the 
ik 7^ 0. Moreover, when one wants to compute even a very simple trace : for example tq,n 



2 



by using Eq. (0), one obtains, [|] steps later, a n — [|] degree complicated trace which 
does not contain quadratic terms any more. Thus, the recursion stops there. In fact, this 
problem can be overcome by a very simple idea : one uses the invariance of traces by cyclic 
permutations to get rid of the n + 1 degree term in Eq. (|^). Then, one obtains relations 
between general traces, and it is thus possible to compute the expectation values of any 
trace in function of the first ones. 

Let us see now how this idea applies to the computation of the resolvent. We denote : 

uq = LO is the usual resolvent. Using the change of variables : 

+ e^— (9) 
z — <P 

we obtain the equation : 

z {U + gz) 00 - U - gz - g - uj'^ + 2cuJi_i = (10) 

Similarly, 

$ ^ $ + e5i<l>5_i^^ (11) 

yields : 

z {U + gz) oJi^i — {U + gz) tq> — g ti-i<i> — u ui^i + cui o-i + cu;_i_i_i = (12) 

and, by the means of similar changes in variables, we have the equations : 

z{U + gz) - {U + 5(2;) ti_i$ - + ct^_i 0-1-1 + ccl^i - cl^ u;_i_i_i = 

(13) 

z{U+gz) uji i-i-i-(t/+fi'2)t-i-i-i$-fi'ti i-i-i$+cu;i i-i-i+cu;_i_i = 

(14) 

These equations alone are not sufficient to compute uj{z). Indeed, if we intend to 
calculate a; (2;), we generate the function u;i_i(z) (Eq. (pUj)). Then, in turn, we generate 
the function co'_i_i_i(z) (Eq. (0)) and so on. 

As for oj functions containing a (i.e. a $^ term) such as ui 0-1, they are easy to deal 
with : we know how to compute traces containing cuio-i = 3^tr5i$^(5_i^3^, will be 
seen as ^tr$^(5_i —^5^ . Then the change in variables : 

AN ^ z — <I> ^ ^ 

$ $ + e5_i^^5i (15) 
2; — $ 

yields (tui-i = tu-i 1 for symmetry reasons) 

g oJiQ-^i + {U + c) oji^i + czuj — c = Q (16) 
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and similar changes in variables lead to the equations : 

guj-i 0-1-1 + U Ci;_i_i_i + aui o-i + czlui^i — cti^, = (17) 

9^1 1-1-1 + UuJi + czu_i_i_i - cti_i$ + cu;_i o-i-i - t$u; = (18) 

But, to compute uj^i-i i-i-i, as mentionned in the comments to Eq. (|^), we have to 
substract two different changes in variables and use cyclicity of traces : 

$ $ + e[$(5_i$5_i(;z - $)-^(5_i$ - (5_i$(5_i(z - ^yH.i^'^] (19) 

yields : 

c {uj-i-i 1-1-1 + uj-i 1-1-1-1 - w-i-i 0-1 - cl^-i 1 i-i) - ^-1-1-1 + tUi-1 = (20) 
This equation, as we know how to compute tui-i, cij_i_i_i, cj-i-i o-i and cj-i o i i-i, relates 

f-i^-l-l 1-1-1 to UJ_i i_i„i_i. 

$ ^ $ + e($5-i$5-i$5-i(2 - <l>)~i$ - (5-i<l>(5-i<l>(5-i(^ - $)"i<l>2) 

allows us to relate similarly uj^i i-i-i-i to tui-i-i-i-i. Then 

$ ^ $ + t{^5^i{z - $)-^(5i<l>(5-i<l>(5i - 5^i{z - $)-i5i$5-i$5i<l>) (21) 

allows us to relate cui-i-i-i-i to dJi-im, and we have cji-i i i i = oj-i i-i-i-i as the roles 
of 5i and 5-1 are completely symmetric. 
Finally, as a result of these operations, we have : 

tu-i 1-1-1-1 + io-i 1-1-1-1 = K{z) (22) 

where K{z) only contains easy to compute uo functions. We can then write an equation 
for uj^i 1-1-1-1 and thus for uo^i^i i-i-i which only involves uj functions that we either 
already know or are able to compute similarly as was done during the two first steps of the 
procedure. 

That way, our set of equations is closed, and we obtain a degree five algebraic equation 
for uj{z). This expression applies to general expressions of U and c. For the exact expression 
of this equation see Appendix A. The equation only contains four unknown parameters : 

U = ^(tr$i), ti-i$ = ^(tr$i$2), till* = ^(tr$i$2$3), h i-i-i$ = ^(tr<l>i$2$i$3) 

(23) 

These parameters are also those that would be involved if we used the renormalisation 
group method [jlO|, 0, H to compute the Potts-3 model. The renormalization group flows 
would relate the conventional Potts-3 to the Potts-3 + branching interactions model, with 
arbitrary U and c; but the presence of tiii$ shows us it would also be related to the dilute 
Potts-3 model, where one has a -i:tr($i + $2 + ^3)^ term. Finally, the ti i-i-i<i> term shows 
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us it may also be related to more complicated quartic models. 

We are now going to derive from our equation the critical behaviour and critical line 
of the model when U = 1 + hti^'^/G and c is a constant. This is the most common type 
of extension of a matrix model to branching interactions. The values of the unknown 
parameters given in Eq. ( 123|) are fixed by the physical constraint that the resolvent has 
only one physical cut which corresponds to the support of the eigenvalues of $. Then, one 
can study the critical behaviour of the model. 

It is easy to look for the Potts critical line. Indeed, the scaling behaviour of the resolvent 
is then, if we denote the physical cut of uj as [a, b] : 

uj{z) ~ (2; — a) 2 when ^ ~ a and uj{z) ~ (2; — 6)5 when 2; ~ 6 (24) 

The corresponding exponent 7^ is — |, which corresponds to the C* = | central charge of 
the model. 

Rather than looking for the resolvent for any values of the coupling constants, it is 
easier to search for the resolvent only on this critical line where the presence of the | 
exponent leads to simple conditions on the derivatives of the algebraic equation. 

We obtain : 

105c' + Ag^ = 

2480625 (-1 - 4c + 43c^) + 296100 c (15 + 113c) h - 692968 h"^ = ^ ^ 

Let us note here that, when h = (no branching interactions) we recover the Potts-3 
bicritical point which agrees with Daul's result : 

_ 2-V47 _ VToE f -3 + V47 \ ^ 
^~ 43 ' ^ ~ 2 V 41 - V47 / 
Thus, we have shown that the resolvent for the model of Potts-3 plus branching interac- 
tions obeys a degree five algebraic equation. We have found the critical line and exponent 
of this extended model. This extends the results of Daul who had only derived the 
position of the critical point and exponent of the conventional model. Finally, let us recall 
that, in a recent paper 0, P. Zinn- Justin obtains independently algebraic equations for 
similar problems. His method, though, does not involve loop equations, and is rather in 
the spirit of ^ . Moreover, it does not address the problem of branching interactions and 
thus overlaps our results only in the case of the conventional Potts-3 model. 



3 The Potts-00 model : 

We are now going to briefly derive the solution for the Potts-00 model, from the equations 
of motion point of view. The purpose of this part is mainly to show the efficiency of our 
method on this c = 00 model. This model was previously studied by Wexler in |Tl| . 
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/ $1 \ 



Let us denote $ = 

V %J 
We shall define the Potts-q partition function as 



and X= (5i±-±5^0V,. 



2 i-^v2 



Z = J rf$e-^^(*) where = 9^ + + c-f^ (26) 

is of order q when g ^ oo. First, let us use the equations of motion to relate 

Let us also denote : 

d(x,y) = —{tT -) (28) 

$ ^ $ + e — — ^—T yields 

x-^y -X 

{x{gx + U)+cy- a{x) - ^) d{x, y) + g - gy b{y) - c a{x) - b{y) {gx + U) = (29) 

Q 

We can get rid of d{x,y) since, when x (gx + U) + cy — a{x) — ^ = 0, d(x,y) remains 
finite, thus g — gyb{y) — ca{x) — b{y) (gx + U) = 0. This is sufficient to relate a{x) to b{y). 
Moreover, the value of b{y) is easy to compute when q = oo. 

Let us briefly summarize this computation : we calculate the value of (trX") in the 
g — > oo limit. 
First : 

(trX") = (tr<l>i...$„) + 0(i) (30) 

q 

(recall that all the play the same role). 

If we now separate the first n matrices from the remaining q — n (with q ^ n), and 
suppose there is a saddle point for the eigenvalues of ^n+iljid^^ then this saddle point is 
(in the g — > oo limit) independent from the matrices $1 , . . . , Then, in this limit, up 
to a change in variables : $fc = U^kU~^, we have n independent matrice $1 . . . Each 
of them has the partition function 

= J e-^(ft^*'.+Ttr*i+ctr*,Ao) (31) 

As tr$i . . . $„ = tr<l>i . . . <l>n, we have (trX") = tr(<li)A^ . . . {^n)Ac where (. . .)ac is the 
expectation value obtained with the partition function (cf Eq. (0) ). 
The matrices $fc , /c = 1, . . . n all play the same part, and (trX") = trA^, thus 

trAJ = tr(<li)" (32) 
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This must give us Ac, provided we calculate ($i)ac iii function of Ac- This is a solvable 
problem, but it is much faster to note that 

Ac = U ^qNxqN 

(33) 

is solution. Thus, ^(trX") = (t$)", and b{y) is simply {y — t^)~^. 

This gives us immediately the solution for a{x) : it obeys a second order equation and 
reads : 

a(x) = ^{x{U + gx) + cU- ^ {x{U + gx) + cUY - A{U + gx + gU)) (34) 

The Potts-oo plus branched polymers model is thus very similar to an ordinary pure 
gravity model. As previously, we compute the parameter t$ by imposing that the resolvent 

3 

a(x) has only one physical cut. The model is critical when a{x) behaves as (x — Xo)2, xq 



being a constant, and the critical point verifies (as in |jTl 



1,1 

9c = ^ and = -- (35) 

Let us note finally that the loop equations method used here is appropriate for the 
renormalization group method of |l^, 0, H . 



4 Conclusion 

In this letter, we have shown that it is possible to solve the Potts-3 and Potts-oo models 
on two-dimensional random lattices through the method of the equations of motion. We 
have obtained a closed set of loop equations for the Potts-3 model, which was thought to 
be impossible. We have shown that the Potts-3 resolvent obeys an order five equation, 
and this new knowledge opens the door to the calculation of expectation values of the 
operators of the model. We have extended the Potts-3 conventional model to Potts-3 plus 
branching interactions, and given the general algebraic equation and the Potts critical line 
of this model. Finally, we have shown our method also applies successfully to another 
Potts model : the Potts-oo model. We hope to generalize soon our method to more general 
Potts-g models, in particular for large-g Potts + branching interactions models. 



A The equation for the Potts-3 resolvent : 

Here is the degree five equation for the resolvent of this model, where W(x) is related to 
u!{x) by W(x) = u!{x) — gx"^ — Ux. 

-24 + 4 _ ^5 g2 _ 12 ti,,^ c^g^-A c^g^ + S tn-i-i^ / + 68 ^ + 

2c^gH^ + 3 g^tl + GOc^U + 2c^g^U - 52 g"^ U + 20 tni<p g'^ U - 20 gt^U - 
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Ac^gH^U -36 U^-3c'^ g^U^ + 36 h^i^ - 36 c"^ g -12 c^U^ + 36 c^gt^ + 

12c^ + 28c^ g X + 2c^ g^ X - 12 g^ x + 8 tm^ g^ x - 36 c'g^t^x-2 gU^x- 

52c^ gU X - Ac^ g^U X + 36ti-i^c^ g^U X - 18c'^ g'^t4,U X + 2c'^ gU^ X + 5Ac^ g^t^U'^ X + 
22c^ gU^ X-2A c^g^x'^-c^g^ a;V8ti_i<^ g^ x^ + 2c'^ gH^x^ + 24 U x^ + 2{]c^ g'^U x^ + 
26 g^t^U x^- 48 c^U^x"^ + 12 g'^U^ x^ + 2Ae'U^ x^ + 2Ac' gx^ + 6c^ g^ x^+A g^ x^ - 
6Ac^ gU x^ + 2c^ g^U x^ + 44 c> gU^x^- 16 g^ x^ + 2Ac^ g^U x^ + A c> g^x^ + {2Ac^g- 
12t,_,^c'g^ + 8tni^c'g^-36c''gH^-2cgH^-A0c''gU-2cg^U + 36t,_,^c^g'U- 
18 g^ U-10 g U^+54: g^ U'^+26 g U^+2A s-24 g"^ x-2 c g^ x+16 g* x+ 

Ac^ g^ t^x-36c^Ux + Ac^ g^Ux + 52c^ gH-(f)U x + 12c^ x + 36 g^U^ x-12c^ x + 
12 U^x-Ac^ gx^ + U c^g^x^ + Uc"^ g^t^ x^ -36c^ gU x^ + 10c^ g^ U x^ + 30 g x^ + 
22 gU^x^- 40 g^ x^ + 60 g^ U x^ + 12 g^ x^ + 18 g^x^ + 2 g^Ux^) W(x) + 
( 12 c^-e c g*+2 g^ t^-12 c' U-A g"^ U+26 c g^ U-12 V^+18 c g"^ U^+ 

12 U^-AA c'gx-2c^ g^ x + 12 cg^t^x + 30 gU x + 20cg^U x + 26c^ gU^ x + 6c^ g^ x'^ + 
2cg^x'^ + 12c^g^Ux^ + 39c'^g^U^x^ + 20c^g^x^ + lAc'^g^Ux^ + c^g^x^ ) W{xf + { - 
12c^ g -2cg^ + AgU^- lOc^ gU + Ag^U + 260^ gU^ + 20 c^g^x + Ag'^x + 12c'^g^Ux + 
18cg'^U'^x + 22cg^Ux^ + Acg'^x^) W(x)^ + { - {c^ g^) + 18cg'^ U + Acg^ x + Ag^ U x + 
Ag^x"^) W{xy + Ag^ W{xf = 

Note that U and c may depend, in the most general case, on ti-i^ and t$2, the latter 
being related to t$ through the equation of motion : gt^2 + (U + 2c)tq> = 0. In this 
article, we have computed explicitely the critical line for the particular case of c constant 
and U = 1 + 1 1$2 . 
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